Quantum trajectories of interacting pseudo-spin-networks 
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We consider quantum trajectories of composite systems as generated by the stochastic unrav- 
eling of the respective Lindblad-master-equation. Their classical limit is taken to correspond to 
local jumps between orthogonal states. Based on statistical distributions of jump- and inter-jump- 
distances we are able to quantify the non-classicality of quantum trajectories. To account for the 
operational effect of entanglement we introduce the novel concept of "co-jumps" . 



I. INTRODUCTION 



In contrast to the ensemble description based on master equations, quantum trajectories make available further 
information about details which are lost by ensemble averaging. The study of quantum trajectories of open systems 
|jj H has therefore found broad application in various fields to analyze fundamental processes. 

It is important to distinguish the quantum trajectories in Hilbert space from the Fcynman paths whereas the 
latter are defined in the underlying classical phase space and contribute with a complex probability amplitude to the 
path integral, the Hilbert space quantum trajectories contribute with a real and positive probability to the ensemble 
■ density matrix Q . 

t-H \ Hence in the Feynman path integral formulation the classical limit can be obtained in a direct way when the interfer- 
ences between the paths basically reduce to the classical path. In contrast the Hilbert space does not have a direct 
classical analogue. Here, we define a Hilbert space-trajectory as "classical" , if it is constrained to orthogonal states. 
The jumps in the Hilbert space imply jumps in observable space like energy or angular momentum: The orthogonal 
states can then be interpreted as the eigenstates of some observable; its spectrum will here taken to be discrete and 
Qh! finite ( "telegraph signal" ) . 

4L ■ While the non-classicality of states has attracted much interest recently J^-§] , the non-classicality of the dynamical 
evolution, in particular on the level of trajectories, has received little if any attention so far. It is tempting to expect 
that an increasing dissipative interaction with the classical environment should make the trajectories more and more 
"classical". One may wonder, however, to what extent this expectation can be verified quantitatively. One may also 
wonder whether there is a relationship between non-classicality of states and of trajectories. 

In this paper we introduce statistical distribution functions to account for the non-classicality of quantum trajec- 
^> , tories, show their properties and their relation to the non-classicality of states. 

This paper is organized as follows: In section |lj we define a convenient operator set for the description of quantum 



networks, the states of which will be discussed in section III. The unraveling of the Linblad master equation into 



single quantum trajectories is summarized in section IV A. In section IV B we d iscus s som e properties of the measure 
of state distance, which are applied to the quantum trajectories in section IV and IV D. The concept of jump- and 



co-jump-distances is exemplified in section [v| for special two-, three- and four-particle states. In section VI we present 
numerical results exploiting this jump-concept for simulated trajectories. We conclude with a brief summary. 

II. CLUSTER-OPERATORS 

We consider a network consisting of N subsystems of n states each. The local states are \p(fi)), p = 1,2, ... ,n, 
/j = 1,2, . . . , N, allowing to introduce the transition-operators P pq (n) = \p(fi))(q(fi)\. These can be combined to give 
the n 2 generators of the SU(n)-algebra, which read for n — 2: 

Ai(M)=A 2 (/i)+iW) (1) 
A 2 (M)=i(A2(M)-P 2 iG«)) (2) 

A 3 (M)=iW)-Ai(M) (3) 
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Ill QUANTUM STATES 2 

Ao(/i)=Pii(M)+A2(M) = l(M) (4) 

They constitute a complete, orthogonal set of n 2 — 1 traceless operators. A corresponding set for the total network is 
then given by the n 2N product-operators: Here we will restrict ourselves to N — 4, n — 2, in which case we have pQl 

Qmikj = A m (4) ® Aj(3) (8 A fe (2) ® A^l) (5) 

The number c of indices unequal zero is the number of subsystems this operator acts on. There are n c — (J^)(n 2 — l) c 
such c-cluster-operators, with < c < N. Operators Q acting on different subsystems commute. 



III. QUANTUM STATES 
A. Correlations 

Any network-operator in a given Liouville space {N, n} can be expressed in terms of such cluster operators: In 
particular, for the density-operator p one finds (again for N = 4, n = 2, cf. Ref. JTo[] ) 

P = 2^ KmlkjQmlkj (6) 

jklra 

with K m lkj = tr jpQmifcj}- 

These expectation- values K = {K m ikj} uniquely specify the state; they decompose into n c c-point correlation func- 
tions. The only c=0-term is -Koooo = tr{/5} = 1. The c=l-terms are the (local) Bloch- vectors, -Koooj = A^ , 

Kooko = A^ etc., which can be found from local (ensemble-)measurements. The c > 1-terms are typically inferred 
from coincidence-measurements (ensemble measurements). 



B. Covariances 



One easily convinces oneself that these correlation functions K m ikj factorize if and only if the state p exhibits some 
product form: For example, if p(4, 3, 2, 1) = p(4, 2) eg) p(3, 1), then 

K m ikj — K m0 ko ■ Koioj (7) 

or if p(4, 3, 2, 1) = p(4) p(3, 2) <g> p(l), then 

K m ikj — K m0 oo ■ K iko ' K n0 ()j etc. (8) 

Subsystems and groups of subsystems which do not factor are called "entangled" (p of the total system is taken to be 
pure and to describe a single network). Without any entanglement, all c-point correlation functions are thus reducible 
to local expectation values (i.e. of type c = 1). 

It is therefore convenient to introduce state-parameters, which describe deviations from this factorization property. 
For this purpose we introduce a supplementary set of cluster-operators 

mlkj 

= AA m (4) ® AA;(3) ® AA fc (2) $ AA^l) (9) 

based on the local "deviation-operators" 

AA ro ( M ) = { X ?{^ - X ™ 1 ^ for m ^ (10) 
\ for m = y ' 

The respective expectation values ( "quantum-covariances" ) 

M m ikj = tr ^AQmlkjp\ (11) 
then also come in different c-classes: For c = 0, Moooo = -Koooo = 1> f° r c = 1> -^oooj = etc., for c = 2, 
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Mookj = Kook] — Kooko ■ K 000 j etc. (12) 

for c = 3, 

Moikj = Kaikj — Kaika ■ X^ — Koioj ■ X^ 

-Kookj ■ A| 3) + 2A, (3) • \™ ■ Xf ) etc. (13) 

The set of expectation values {X^\ M m jjy} can alternatively be used to specify the network-state. With all X h ^ = 
we obviously have M m ikj — Kmikj- The factoring properties of K m ikj carry over to M m ikj- In particular, under the 
condition as for eq. (||) we get M m ikj = M m ooo ■ Moiko ■ Moooj = 0. In general, any specific M m ikj is zero, if at least 
one individual subsystem entering with a local operator-index ^ factors out. 

C. Entanglement measures 

For a product state as of eq. (0) we have M m ikj — M m ofeo ■ Moioj ^ 0. By substracting all possible partitions 
(c = c\ + c 2 + . . ., Cj > 2) we introduce (here for c = 4 = 2 + 2): 

M m ikj = M m ikj — M m i()oMo()kj 

— M m okoMoioj 

- MmOOjMoiko (14) 

which is thus zero for any product state. For c < 4 we obviously get M = M; for c > 4 this connection scheme is 
easily generalized. 

To quantify entanglement on the total network-level (N = 4) we use (cf. Rcf. fnj) 

3 2 

/3(4,3,2,1)= (Mmikj) (15) 

as well as corresponding sub-space measures like 

3 2 

/3(2,1)= ^ (Mooy) (16) 

Note that (3(2, 1)^0 indicates any entanglement M 00k j ^ between subsystems (2) and (1) only: /3(4, 3, 2, 1) 
could still be zero. These /3-functions can be used instead of the "entropies of entanglement" (=entropy of the 
respective reduced density operators) ||[|; the M-terms are easier to calculate and, furthermore, can be made basis 
of approximation schemes (see below). 

IV. QUANTUM TRAJECTORIES 
A. Stochastic unraveling 

The evolution of open quantum systems is usually approximated by the Lindblad Master equation, which is Marko- 
vian. Here, the influence of the environment is specified by so called environment operators L s and the corresponding 
damping rates W s : 



-p + l -[H, p] = ^(--W s {LtL s p + P LtL s } 



dt' h l ,ri ^ y 2 

S 

+W s L s pLf) (17) 

Those operators L s play a crucial role in the stochastic unraveling of the master equation. The coupling to the 
environment leads to individual quantum jumps between which there is a non-unitary continuous evolution [jjj. The 
jumps are generated by the last term of the right-hand side whereas the first two terms can be combined with the 
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Hamilton operator into a non-Hermitian effective Hamiltonian responsible for the continuous inter-jump-evolution. 
The probability for a quantum jump of type s after the time interval St is given by 

p s = W s tr{L s pL+}5t (18) 

We suppose that these projections can be expressed in terms of, in general, non-Hermitian operators, L s = P s (p) 
which are taken to act locally on one of the subsystems (p) , 

P s (p)pP s (p)+ 

P = — rz : T ( 19 ) 

tr|P s (/i)pP s (/i)+} 

For the characterization of these pure-state-trajectories the timing of jumps plays a central role because it may give rise 
to measurable events and to count- and waiting-time-statistics jyj. Unfortunately, however, these give only indirect 
evidence for the non-classicality of trajectories (like, e.g., ant i- bunching ) . We therefore propose to supplement the 
analysis by directly referring to the motion in Liouville-space. 

The characterization of single quantum trajectories with respect to classicality measures could be done in different 
ways. In addition to the possibilities presented below one may think of calculating the distribution functions of local 



as well as non-local coherence measures (cf. Sect. (QIC)). However, it should be noted that these properties depend 
on the basic operators chosen for the state description. In the case that the coupling to the environment leads to the 
built-up of states which do not happen to coincide with eigenstates of the local operators A3 (p), this approach will 
not show a classical limit (i.e. a and/or (3 remain unequal zero). 



B. Measures of state distance 

There have been different proposals for defining a metric for density matrices (see, e.g. the Bures metric [fl2|j , fidelity 
[ 13 1 or mutual information H). For the non-orthogonal Glauber-states |a), e.g., a "distance" d has been proposed 
[14| with Kale/)! 2 = exp{— cr}. In this paper we use a measure, D, for the distance between two arbitrary (generally 
mixed) states p and p' according to 

D%,=tv{{p-p'f) (20) 

which is, independent of the dimension of the Liouville space, between and 2. The maximum (squared) distance of 
2 applies to orthogonal states. In the case of pure states (p = |^ r )(^|) D 2 ~~, can be rewritten as : 

^,=2(1-|W)| 2 ). (21) 

It is easy to show, that Dpp* satisfies the metric properties |l5| in the Liouville space C, i.e. 

Dpp' > for all p, p' of £, 
Dpp' =0 if and only if p = p' , 
Dpp' = Dp, p for all p, p' of £, 

Dpp' < Dpp" + Dp"p' (triangle inequality). (22) 

This measure D 2 ~~, can directly be expressed in terms of the SU(2)-parameters, K m u~j, namely as the squared length 
of the difference vector between K= {K m ikj} and K'= {K' m ikj}i 

1 3 2 
Djp' = ^2 { K mlkj - K'mlkj) ( 23 ) 

j/,fc,/,m— 

This concept of state distance can easily be generalized to distances defined on reduced state spaces: Observing that, 
e-g-, 

tr {4,3,2} = Aj(l)2 5 m0 5i 5ko (24) 

(here tr{ M } means trace-operation within /z-subspace only), we find for the reduced density operator of subsystem (1), 
say, 
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= tr {4 32}/5 = ^^000,^(1) (25) 

3 

and 

K!p'}) 2 = 1E(^oo,-^oo,) 2 (26) 
Correspondingly, the distance as seen from the subsystems (4, 3, 2) is 

Hs?) 2 4^ (K ™' fco - ^*° )2 (27) 

k,l,m 

Then we have by inspection the inequality 

C. State-distance distributions 

In order to characterize quantum trajectories we introduce various types of state-distances: the "jump distance", 
by inserting into eq. (|2l]) the state right before and after the jump (in analogy to the jump distance of the Brownian 
motion in classical physics the "inter-jump-distance" as the distance between the final state of the last jump 

and the initial state of the following jump, and the state distance for a specified time interval r during the evolution 
of a given quantum trajectory, 

Dl = tr{{p{t)-p{t + r)f}. (29) 
Finally, we will be interested not only in the jump distance of the total system, but also of parts of the system in their 
respective reduced space p. Therefore we use ^D^l,jJ as defined in eq. ( |26| , |2"7| ) where p (p') is the total density 

operator before (after) jump. In this way we can test to what extent a projection in subspace (z/; v ^ p) affects the 
reduced state of subsystem p ("co-jumps"). 

Sampling over one individual trajectory we find the corresponding distribution functions, f(D 2 ). These are nor- 
malized: 



[ f(D 2 )d(D) 2 = l. 
Jo 



(30) 



D. Co-jumps and entanglement 

One easily shows that for any complete (POVM)-type measurement jnj, the ensemble-averaged co-jump must be 
zero: For this purpose we write, for the (ensemble) density-operator after measurement in (2,1), e.g. 

/ }' = ^A(2,1) / 5P S + (2,1) (31) 

5 

with J2 S -P,s + (2: 1)-Ps(2, 1) = 1 and consider (overlines indicate ensemble-averaging) 

K mlQ0 = tr {/5'Qm/Oo} (32) 



tr|^](p s (2,l)pP+(2,l)Q m . 
^^2(P+(2,l)Ps(2,l)) pQ 



(33) 



(34) 



tr 



{pQmZOoj = -KmiOO (35) 
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Here we have made use of the fact that P s + (2, 1) and Q m wo commute as they act on different sub-spaces. We conclude 
that 



( D w\Y = 4 E ( R -'oo - <.,oo) 2 = (36) 
il 



This means that ensemble quantum mechanics is local in an operational sense: Measuring in some subspace, here (2,1), 
does not have any influence outside this subspace. An analogue statement holds for respective unitary transformations. 

However, as is known since the famous EPR-experiments jTJj individual measurements leading to new information 
may violate this locality. To show this we consider 

p' = -Ps (2, 1) / 5P+(2,1) (37) 

Ps 

where p s — tr |p+(2, 1)P S (2, l)p|. Based on the same arguments as before we obtain 

K'mlOO — KmlOO H trj p ( QmlOO — 1-KmiOo) 

xP+(2,l)A(2,l)} (38) 
Now let P+P s = \ (1(1) + A 3 (l)) = P 22 (l): Then we get 

K'mlOO ~ KmlOO = 7) — {Kml03 — KmlOQ ' ^0003) (39) 
^Ps 

The right hand side is zero if -RT m ;o3 = K m ioo ■ ^0003, i-e. if subsystem (1), which is measured, has no entanglement 
with the subsystem (4,3). Otherwise, 



{ D ipP}) 2 = h E ^™' 00 ~ <<oo) 2 ^ (40) 



2 1 

2 

which is thus non-zero also for the ensemble: Comparing with eq. (|36j) we note that 



(K mWQ - K' mm ) 2 ^ (K ml00 - K' rnm ) 2 . (41) 



V. CO-JUMP PROPERTIES OF MODEL STATES 



Our intention is to use the concept of co-jumps for the characterization of single quantum trajectories. As these 
trajectories always connect pure states, "mixed states" only appear for reduced subspaces. Any non-zero entropy of 
such reduced density matrices is due to entanglement and not due to our incomplete knowledge. Co-jumps will occur 
also within reduced spaces; such situations will be included below in a formal way. 

In the following we perform individual "measurements" based on the specific operators (cf. eq. (|H 

A(M) = \ (l(M) +Ai(/x)) ,» = 0,1,2,3 . (42) 

Note that these measurements correspond to only one outcome each. In the case of coincidence measurements we 
take the dyadic product of such single particle operators. According to eq. (23) and eq. ( |26]j27j ) we then calculate the 
jump- and co-jump-distance. 

A. Two-particle state 

We first consider the completely mixed state, the EPR-state ("cat-state") |EPR) = -j= (|12) — |21)), and the 
"Werner-state" B 
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p x = (l-:r)ii + a;|EPR){EPR| ,0<a;<l. (43) 



The maximum entanglement for the latter is reached for x = 1 (EPR-state, /3(2, 1) = 3). Table |V A| shows results for 
x = and x = 1, conditioned by the respective projection, i. Jump and co-jump are independent of i, confirming the 
"isotropy" of the EPR-state and of the mixed state. 

In Fig. [j] one can see the transition from the mixed state to the EPR-state. In contrast to D, the co-jump is 
linear in the whole region from x — to x = 1 : The state (^3|) reacts in a local way for x — only. 

The non-locality, to be sure, could still be "explained" by a local hidden- variable theory || unless x > 1/V% 
(violation of Bell- inequalities) or x > 1/3 (violation of separability condition 0), respectively. These interpretations 
have to assume, though, that mixed states can be treated as classical mixtures (i.e. resulting from our ignorance 
|jl9f rather than from undefined properties due to entanglement with other subsystems). Stochastic modeling deals 
with single networks and individual subsystems (additional knowledge!); co-jumps within the latter then reflect their 
("objective") change as induced by a distant observation. 

TABLE I. (Squared) jump distances for the specific iV=2-particle state p x (see eq. (|43[)) (projection by P^) 

i (D^) 2 (Df 

1 1 0.5 1.0 

1 2 0.5 1.0 

1 3 0.5 1.0 

1 0.0 0.25 

2 0.0 0.25 

3 0.0 0.25 
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B. Three-particle state 



As a representative for the 7V=3-case we consider the cat-state (GHZ-statc, /3(3, 2, 1) — 4, (3(2, 1) = 1) 

|GHZ) = -±= (|111) + |222» 



(44) 



and study projections on particle 3, and on particles 2 and 3 in coincidence. In contrast to the EPR-state, the GHZ- 
state is not rotationally invariant, the t=3-direction play s a special role which can be seen by comparing the result 
for = (0,1) and (0,3) or (1,1), (3,3) (see tabic |VB|). 



TABLE II. (Squared) jump distances for the iV=3-cat-state (projection by P\ 
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One should note that the reduced state p(3, 1) could be written as a mixed product state (i.e. separable in the sense 
of Ref. fif]), which, nevertheless, leads to a co-jump in (1) induced by subsyste m (3 ). Only without further knowledge 
could this effect be interpreted as being due to classical correlations, cf. Sect. V A[ 



C. Four- particle state 



The special four-particle "cat-state" ((3(4, 3, 2, 1) = 12, (3(3, 2, 1) = 0, (3(2, 1) = 1) 



I* cat) =-J=(|llll) + |2222)) 



(45) 



shows a co-jump behavior which is, contrary to the 7V=3-cat stat e, di fferent also for (i,j) = (1, 1) and (1,3). The 



pertinent co-jump properties of this state are summarized in table V C 



TABLE III. (Squared) jump distances for the A r =4-cat-state (projection by f^ (3) ® P w ) 
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VI. JUMP STATISTICS OF SIMULATED TRAJECTORIES 
A. Hamilton-model 

In our simulations we consider an open quantum network consisting of pseudo spins (two-level systems) which 
interact with each other, with external electro-magnetic fields and the environment. In rotating wave approximation 
the local Hamiltonian of node /i can be expressed as (fi = 1, 2, 3, 4) 

H(p) = \ <f" A 3 (/x) + \ 9* ■ (46) 

The external fields are characterized by the Rabi-frequencies and the energy difference between the laser photons 
and the spin energy is denoted by S^. 

The interaction between the spins will be given by the non-resonant coupling, 

H(m,v) = -\cf v) A 3 (m) ® AsM , (47) 

which leads to an energy shift of one subsystem depending on the state of the other. In the following we will call 
a network "homogeneous" , if all model parameters (including those for the bath coupling) are invariant under any 
subsystem-index permutation. 

B. Bath model 

In addition to the Hamiltonian model, we have to specify the coupling to the environment leading to decoherence 
effects. For simplicity we will take into account only dissipation of energy into the bath with rate Wj* and operator 

— Pi2{n) an< 4 for the case of a bath with finite temperature also excitations out of the bath with rate and 

operator Lf(fi) = Pai(fi) 

In the case of a bath temperature T = there are only jumps into the ground state, from where only the coherent 
laser field can drive the state out again. Therefore the distribution functions for jump- and inter-jump distances are 
virtually identical. 




FIG. 2. Distribution function f(D y 9 t) for the jump distance D of a homogeneous network TV in interaction with a bath at 
T — (damping rate = W , — g — coupling parameter to the coherent driving field, coupling Cj£ = 20, detuning 
5" = -20 for TV = 2 (-40 for TV = 4) ) 

For small environmental influence the broad distribution reaches its maximum at = 2, due to the fact, that 

jumps (photon-emissions) most likely start from the upper level in which case the jump into the ground level has 
(squared) distance 2. Increasing the ratio W/g raises the probability for jumps already for incomplete excitation, 
which results in a (squared) jump distance smaller than 2. In the high damping limit, the jump distance peaks at 
zero and the distribution becomes narrow (see Fig. 0). 
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The jump-distribution function f(D 2 ) for homogeneous networks shows a scaling behavior which is only weakly 
dependent on size N (Fig. |J). Contrary to naive expectation, the limit N — > oo thus does not necessarily mean a 
classical limit. For W/g 3> 1 the dynamics is virtually suppressed. 

The distribution changes completely, if the bath also causes excitations (bath temperature T ^ 0): Now, classical 
trajectories result, if the jump distance distribution peaks at D 2 = 2 and becomes very narrow, while the inter- 
jump-distances tend to zero (Fig. || and Fig. |J). This happens for W/g 3> 1; it will always happen for incoherent 
driving. 




FIG. 3. Distribution function f{Dy 9 ,) for the jump distance D for a network N — 2 in interaction with a bath of high 
temperature (damping rate = W!/ = W, t/ M = g = coupling parameter to the coherent driving field, renormalisation 
coupling C R = 20, detuning <5 M = -20) 




FIG. 4. Distribution function f{D%^) for the inter jump distance, D, for a network N = 2 in interaction with a bath of 
high temperature (parameters see Fig. 



C. Influence of entanglement 



We may truncate the expansion of correlation functions at a certain level, c, i.e. we neglect all orde rs of e ntanglement 
M m ikj above this specified level. This leads to a reduction of relevant state parameters (cf. Sect. IIIB), but at the 
same time to non-linear evolution equations. (For c = 2, we replace, e.g. -f^oifej by the various factors as obtained 
from eq. (13) with Moi/y =0.) As this factorization does not correspond to a concrete partitioning, those equations 
are not guaranteed to remain consistent: To the contrary, pure states (of closed systems) do no longer remain pure, so 
that quantum trajectories based on this strategy lack a definite interpretation as a single system quantum evolution. 
Only in the presence of sufficient damping such deficiencies are negligible. 

Alternatively, we may carry out simulations for which the state is taken to factor into concrete partitions, like, e.g., 
/5(4, 3, 2, 1) = p(4, 3) <8>/5(2, 1). As a consequence, all entanglement terms, M m ikj, between those partitioned subgroups 
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disappear. In this case the equations of motions remain consistent, i.e. the density matrix for a closed system keeps its 
desired properties and a pure state stays pure. The comparison of these simulations with the exact ones then allows 
to isolate the effect of various entanglements M m ikj on the level of individual quantum trajectories. 

The resulting quantum trajectories are characterized as introduced for the exact ones: The respective distribution 
function for the jump distance is plotted in Fig. where all M-terms have been set to zero, a factorization into 
individual spins. 

As we can see, increasing damping rates decrease the difference between the distributions of the exact and this 
factorized simulation, respectively. In the large damping regime the latter can be understood as a good approximation, 
which thus allows for an enormous reduction of the number of relevant state variables. 




FIG. 5. Distribution function /(D^y,) for the jump distance of N pseudo spins (T = 0), neglecting all entanglements 
(parameters see Fig. ^ 

For homogeneous systems adequate factorizations should also be permutation-symmetric with respect to subsystem- 
indices. The only candidates then are truncation schemes. In Fig. |g we show an example for N — 4 where all entangle- 
ment terms beyond second order (c = 2) are set to zero. Nominally, this approach should constitute an improvement 
over the results of Fig. ||, for which all entanglement has been neglected. Qualitatively, the distribution of Fig. ^ is, 
indeed, between the exact and the non-entanglcmcnt-rcsult. 




FIG. 6. Distribution function f{D^ 9 ,) for the jump distance of N = 4 pseudo spins (T = 0), neglecting all entanglements 
beyond second order (parameters see Fig. ^) 

However, its significance is difficult to assess, as for W/g 3> 1 entanglement becomes negligible, anyway, while for 
W/g — > the truncation scheme becomes inconsistent. In this region the jump distance may assume values greater 
than 2, which has no physical meaning. 

On the other hand, there are examples in which the suppression of entanglement can affect quantum trajectories 
even in a qualitative way. This is demonstrated in Fig. ffl for an inhomogeneous N = 2 network. Whereas the exact 
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evolution develops an apparently classical telegraph signal between "light" and "dark" periods ( "Zenon-effect" , ) , 
which appear to be discontinuous only on a large time scale, these transitions vanish, if entanglement is suppressed! 
Fig. [j] shows a section of the quantum trajectories with and without entanglement and the distribution function of 
the state distance D"^., referring to a suitable time scale (r > 1/W). The exact simulation has a peak for D\ — > 0, 
which is missing in the approximation. 
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FIG. 7. Trajectories of a two spin-system and its distribution function f(D T ) for the exact simulation (a) and with neglect 



of entanglement ((b); parameters: Wf = 1.0, g 1 = 0.7, g 2 = 0.03, 8 1 



-20, Cr = 20, t = 1.5 



D. Co-jump distribution 

We finally address the non-locality in quantum trajectories: Using the concept of co-jumps we investigate the 
distribution function of the state changes in the respective non-projected (non measured) part of the system. This 
part changes stochastically from jump to jump; we therefore denote the respective co-jump-distance as D^ red \ Fig. 
H shows this entirely non-local effect for N = 2 and N = 4. As one can see, entanglement is still present even in 
the regime of W/ g ~ 4, where the jump- and inter-jump- distribution functions already indicate a rather "classical" 
behavior. While the classicality of states implies the classicality of trajectories, the inverse is not true, in general. 
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VII. SUMMARY AND CONCLUSIONS 



We have investigated quantum trajectories of networks with up to N = 4 subsystems in terms of various statistical 
state distance distributions. We have found that jump- and inter-jump distributions allow to describe the classical 
limit as well as deviations from this limit in a transparent way. Other measures (like entanglement measures) provide 
supplementary information. 

The classification of quantum trajectories concerning non-classicality has to be distinguished from the study of non- 
classicality of states. Dynamical properties of open quantum systems have been investigated in terms of state changes 
due to quantum jumps and due to the continuous evolution. The co-jump distribution may indicate non-classicality 
of states despite the fairly classical trajectories. 

For this analysis we found it necessary to abandon the "ignorance interpretations" of non-locality, which are often 
applied in quantum-information-scenarios. Non-locality with respect to individual pure-state trajectories is, instead, 
introduced in an operational way, namely in terms of co-jumps. 

In general, large quantum networks do not necessarily behave "classical". In the strong damping limit, "classical 
trajectories" ("telegraph signals") result (for a bath temperature T > 0) which means that the state-space available 
to the quantum network is extremely compressed. This can be exploited for approximation schemes. However, some 
caution should be exercised, as truncation schemes (factorization beyond a certain level c), can lead to inconsistencies 
for small damping. 
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